Landau levels and the Thomas-Fermi structure of rapidly rotating Bose-Einstein 

condensates 
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We show that, within mean-field theory, the density profile of a rapidly rotating harmonically 
trapped Bose-Einstein condensate is of the Thomas-Fermi form as long as the number of vortices 
is much larger than unity. Two forms of the condensate wave function are explored: i) the lowest 
Landau level (LLL) wave function with a regular lattice of vortices multiplied by a slowly vary- 
ing envelope function, which gives rise to components in higher Landau levels; ii) the LLL wave 
function with a nonuniform vortex lattice. From variational calculations we find it most favorable 
energetically to retain the LLL form of the wave function but to allow the vortices to deviate slightly 
from a regular lattice. The predicted distortions of the lattice are small, but in accord with recent 
measurements at lower rates of rotation. 
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How very rapidly rotating Bose-Einstein condensates 
carry angular momentum remains a fundamental ques- 
tion of many-particle physics. Ho 0, noting that the 
Hamiltonian for a rotating gas in a harmonic trap is sim- 
ilar to that for charged particles in a magnetic field, ar- 
gued that for rotational angular velocities just below the 
transverse trap frequency, u>±, all particles would con- 
dense into the lowest Landau level (LLL) of the Coriolis 
force. Motivated by this insight, Schweikhard et al. 
(following earlier work in have recently achieved rota- 
tional angular velocities f2 in excess of 0.99 uj±, at which 
the cloud contains several hundred vortices. The experi- 
ments have reached the "mean-field quantum Hall" fron- 
tier at which TiVL becomes comparable to the interaction 
energy, gn, where n is the particle density, g = Airfi 2 a s /m 
is the two-body interaction strength, m is the particle 
mass, and a s the s-wave scattering length. 

Employing a quantum-Hall like condensate wave func- 
tion with only LLL components, 
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- where the rotation axis is along z, Q = 
(~i are the vortex positions, r = (x,y), d± = y/hj mu>± 
is the transverse oscillator length, and the constant 
normalizes / d 2 t|</>lll| 2 to unity - Ho predicted that for 
a uniform lattice (at each height), the smoothed density 
profile of a trapped cloud would be Gaussian: 
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Here the superscript "u" denotes that the vortex density 
is uniform, 1/er 2 = — Ttn v , n v is the vortex den- 

sity in the plane perpendicular to the rotation axis, and 
(. . .) denotes an average over an area of linear size large 



compared with the vortex separation but small compared 
with a 0. On the other hand, Ref. Q - adopting a more 
general wave function, as in that is a product of a 
slowly varying envelope describing the global structure of 
the cloud and a rapidly varying function describing the 
local properties of individual vortices - found that for 
a uniform array of vortices and for Na s /Z 3> 1, where 
Z is the height in the z direction, the density profile 
is a Thomas-Fermi (TF) inverted parabola. In current 
experiments ; 2f , Na s /Z is of order 10-100, and the trans- 
verse density distributions are in fact better fit with a 
TF profile than a Gaussian. Theoretical evidence for the 
Thomas-Fermi profile has previously been found in nu- 
merical studies 

As first pointed out by A.H. MacDonald (see Ref. 37 of 
0) and stressed to us by N. Read 8|, one can achieve sig- 
nificant changes in the density while maintaining a con- 
densate wave function made only of LLL components, if 
one allows the vortex lattice to relax from exactly tri- 
angular, an effect not considered in Refs. 0, 0|. We 
show here that when the lattice is permitted to relax 
only slightly from uniformity, the TF form is generally 
obtained for Na s /Z > 1 - Q/ujj_. Thus TF behavior 
extends to much lower densities than previously realized. 
This effect comes about through a modification, with in- 
creasing gn, of the condensate wave function, 'F, from 
iV^^LLL- A modification of 'F that allows the density 
distribution to spread out tends to reduce the interaction 
energy: for infinitesimal modification, 'F ~ 
the change in the interaction energy is linear in <5\l/, and, 
with the proper choice of the phase of 8^>, is negative. 
Such modification costs kinetic and trap energy, which 
however is quadratic in for infinitesimal modification. 
Thus the state Q for a uniform array of vortices is always 
unstable towards smoothing out of the density distribu- 
tion. The resulting change in the density profile can be 
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large even though the distortions of the lattice are small. 

We consider the properties of a rotating cloud in a har- 
monic transverse potential V(r) — mw^(i 2 + y 2 )/2. For 
simplicity, we treat mainly the two-dimensional problem, 
and generally set % = 1. The approach adopted in Ref. |j| 
was to express quantities as sums over Wigner-Seitz cells 
for single vortices. This is cumbersome, and in this Letter 
we adopt a different trial condensate wave function which 
is better suited for rotational angular velocities close to 
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The function h, which we assume to be real and slowly 
varying on the scale of the intervortex separation, modi- 
fies the LLL components, as well as admixes higher Lan- 
dau levels. The vortex lattice need not be uniform. The 
wave function @ is much more convenient for calculating 
the kinetic energy than is the ansatz used in Ref. 0. 

The admixture b^^ of the higher Landau level (/i, v) in 
the wave function © , where \i is the angular momentum 
index and v the level index, is of order d\ {d v hjdr v ) r ~ 
(d±/R) v 0(0) - 1) C°, where C° is the amplitude of the 
level [i in </>lll, and r M = d±y/Jt is the peak of the LLL 
wave function n The amplitude for admixture of 

1/2 

higher Landau levels is thus of order d±/R ~ 1/N V rel- 
ative to the modification of the LLL contribution, where 
iV v is the total number of vortices in the cloud. 

The energy per particle of the condensate in the rotat- 
ing frame is E' = E — HL Z , where E is the energy in the 
non-rotating frame and L z is the expectation value of the 
angular momentum per particle about the rotation axis. 
Following Ho Q, we write 



E' = (lua_ - Q)L Z 
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where uj± = uj±z, and $2D is the effective coupling pa- 
rameter in two dimensions. If the system is uniform in 
the z direction, c/2D = Ng/Z, where Z is the axial ex- 
tent of the cloud, while if the system in the z direction 
is in the ground state of a particle in a harmonic po- 
tential of frequency u> z , then g 2 D = Ng/d z V2~Tr, where 
d z ee (Ti/mw z ) 1 ' 2 E3. 

Because the higher Landau levels in the variational 
wave function @ have probability ~ l/N v , the angular 
momentum per particle is 
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plus terms of order 1/N V , which we neglect. Furthermore, 
since h varies slowly in space, we shall use the averaged 
vortex approximation as in to write 



E' ~ Q 
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Here b = (I'/'lllI 2 ) 2 ^ 1-158 describes the 

renormalization of the effective interaction due to the 
rapid density variations on the scale of the vortex sep- 
aration HIH03- 

Uniform lattice. We begin with the ansatz @ for a uni- 
form vortex density. The first term in braces in Eq. JfjJ, 
the extra kinetic energy due to admixture of excited Lan- 
dau levels, scales as 1/-R 2 , as does the interaction energy. 
If <?2D *C ?i 2 /m, the first term suppresses spatial varia- 
tions of h, and h 1, i.e., the density profile is Gaussian. 
This condition is simply that Na s /Z <C 1, where Z is 
the effective size in the z direction, bounded below by 
d z . In the opposite limit, <?2D 3> ft 2 /m, the extra kinetic 
energy is unimportant, and the optimal density profile 
is obtained by minimizing the second and third terms in 
the integrand, which results in a TF density profile, in 
agreement with the considerations of Ref. |j| . 

To obtain illustrative results we employ a variational 
trial function that interpolates between the Gaussian and 
TF (G-TF) forms; exploiting the fact that lim ce _ +00 (l — 
t/a) a — e _t , as in 13j, we take 
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for < r < y/aL, and h = otherwise. Expression 
describes both the Gaussian (a — * oo) and TF (a — 
1) regimes. The number of particles per unit area is 
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thus A 2 = (cr 2 /L 2 )(l + l/a). We refer to 



Eq. for h as the G-TF form. 

Substituting Eq. J7J into |JB}, we obtain 
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Minimization of _Bq_ tf with respect to a and L yields 



L- = <l\ ( 1 - — 

0J\ 
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= L 2 /(l + 2/a) 



(9) 



(10) 



where the dimensionless parameter k ee mbg^D / {2'kTi 2 ) 
determines the strength of interparticle interactions. 

The optimal value of a, determined by minimizing 
E' G _ TF with respect to a, obeys the quartic equation, 
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a 4 - 2(1 + 4A)a 3 - 12Aa 2 + 2(1 - 3A)a - 1 — A = 0, where 
A = k^ 1 . The only real and positive solution is 



2 i 2(1 + 2A) (3277 - 1) 
2a = | 2 + 32r7 + 3-23773 + ' 



y+ + I677- 3 • 23773 



+ (1 + 4A)+ Jl + 1677-3- 2377- 
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where 77 = A(l + A). We note that a is independent of 
O. In the weak (k — > 0) and strong interaction (k — > 00) 
limits, 



a ~ 8/k (« < 1) , 
a ~ l + (3/2i)«;-3 + 



(k ! > 1' 



(12) 
•(13) 



The shape index a, Eq. (|lll) decreases from infinity 
in the absence of interaction to unity as K increases 
from zero to infinity, and the density profile of the cloud 
changes from a Gaussian to an inverted parabola. With- 
out the minimization with respect to a, we can describe 
the cloud by assuming a shape of the density profile cor- 
responding to a given value of the shape index, in terms 
of which L and a are given by Eqs. JHJ and IjlOJI . If we 
let a — > 00, we reproduce a Gaussian density profile as 
in the case without the modulating function, but with 
optimized width. 

We now compare the results for the G-TF profile with 
a optimized and the profile with a — > 00. We write the 
energy for the Gaussian profile as 



n + V2 luj_ (1- — 



(14) 



obtained by taking a — > 00 in Eq. J5J). The relative 
energy difference (E' G _ TF — E' G )/(E' G — SI) between the 
G-TF and Gaussian cases is shown in Fig.^ In the limit 
k — ► 00 , the energy is the TF result, 



E' T¥ ~ fi + -a; j 



si 



K2 . 



(15) 
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and (E' G _ TF - E' G )/(E' G - SI) converges to (2 3 / 2 /3) 
-0.0572, independent of SI (FigJTJ). 

In the experiments in Ref. [2j, where uj z — 2tt x 5.3 
Hz, N — 1.5 x 10 5 , and a s = 5.6 nm for the triplet state 
of 87 Rb, one has Na s /Z > 26 at SI = 0.989 uj ± , typical 
of the highest angular velocities achieved so far. For a 
uniform density in the axial direction, with 52D = gN/Z, 
we find k = 2Na s b/Z > 52, while for a Gaussian density 
profile in z, for which <?2D = gN/d z V2n, we obtain n = 
2Na s b/d z V2n ~ 1.4 x 10 2 . 

In Fig. [21 we plot the normalized density profile v(r) = 
h 2 (\4)LLh\ 2 ) at SI = 0.99 lo ± and k = 100 for the G-TF 
case in addition to v{r) for the two extreme cases of 
the Gaussian form with a — > 00 and the inverted TF 
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FIG. 1: Difference between the energies E' G _ TF calculated in 
G-TF with optimized a, Eq. 1111 . and that for the Gaussian 
case with a — > 00, measured relative to E' G — Q. The curve 
converges to (2 3/2 - 3)/3 ~ -0.0572 as K — > OO. 
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FIG. 2: Global density profile v{r) = /i 2 (|<?!)lll| 2 ) of the 
cloud for Q — 0.99 lu± and n — 100. The dashed line is for 
the optimized a given by Eq. Hill , and the dotted line is the 
Gaussian case. The solid line is the TF profile, Eq. itTTt . 



parabola; see below. The density profile for the G-TF 
modulation (dashed line) is closer to the TF form than 
to a Gaussian. Figure 2] indicates that the relative en- 
ergy reduction compared with E — CI for the Gaussian 
approximation is only ~ 4.95%. 

Distorted lattice. Next we consider an LLL wave 
function with a non-uniform vortex density, and take 
h(r) = 1. Equation © is still valid, and the energy, 
given by 

E' = n + J d 2 r{(^-Sl)^<|0 LLL | 2 ) + ^(|^ LLL | 2 ) 2 }, 

(16) 

is minimized by the TF profile 
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where /i is the chemical potential. Such a solution is 
possible only if one relaxes the constraint of a regular 
lattice; otherwise the solution is a Gaussian, as in pj. 
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The energy is given by the TF result (|15fl . Since the 
term in Eq. <0 involving dh/dr is positive definite, it is 
clear that it is energetically favorable to create deviations 
of the density profile from |<^lllI 2 by deforming the lat- 
tice rather than by exciting higher Landau levels. These 
conclusions have recently been confirmed by numerical 
calculations by Cooper et al. . Excitation of higher 
Landau levels will reduce the energy of a distorted lat- 
tice still further, but it may be shown that this is of order 
(gn) 2 /h£l per particle, which is smaller than the terms 
we have retained. 

Remarkably, the energy and TF profile are valid for all 
interactions strong enough that the averaged vortex ap- 
proximation is valid. This holds provided the size of the 
cloud is large compared with or alternatively that 
the number of vortices in the cloud N v ~ R 2 jd\ is large 
compared with unity. This condition is equivalent to 
Na s /Z » 1 — Q/lj±, i.e., that the energy, U(uj± — fi), 
to excite higher angular momenta [i in the LLL, be small 
compared with the mean field energy in the situation that 
all particles are in the transverse ground state of the trap. 
For fl close to lu± , this condition for a TF profile is much 
weaker than that for a non-rotating cloud. By contrast, 
a TF profile with a uniform lattice requires Na s /Z ^> 1. 

Modifying the amplitudes of lowest Landau level com- 
ponents in the condensate wave function, without intro- 
ducing excited Landau level components, is equivalent 
to changing the positions of the vortices. In the LLL 
wave function , the positions of vortices determine the 
smoothed density distribution generally through 0, 



-V 2 lnn(r) 
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-2- + 7rn v (r). 
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This equation allows us to estimate the displacement of 
the vortices from a triangular lattice for a general density 
distribution. In particular, if we assume a TF distribu- 
tion, n(r) ~ 1 — r 2 /R 2 , where R is the radial extent of 
the cloud, then Eq. (fTH|l implies [l5| 
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The second term is of order \/N v compared with the 
first, since iV v ~ R 2 jd\. The corresponding mean dis- 
placement Sr in the radial direction of the vortices at 
radius r is then 8r/r = (d 2 L /2R 2 )/(l - r 2 /R 2 ) ~ 1/N V . 

Thus in the LLL limit very small distortions of the 
lattice can result in large changes in the density distribu- 
tion. For lower rotation rates, Sheehy and Radzihovsky 
[H| have recently demonstrated that the vortex density 
obeys an equation of the form (|19fl but with a coefficient 
of the second term which depends on the ratio of the vor- 
tex separation to the vortex core size. Such distortions of 
the vortex lattice have been measured experimentally at 
relatively low rotation rates 0, and are in good agree- 



ment with theory [lq|. It would be valuable to extend 
the measurements to rapidly rotating condensates. 
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